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Recently T2K gives hint in favor of large reactor angle S13 . Most of the models, with tri-bimaximal 
mixing at the leading order, can not reproduce such a large mixing angle since they predict typically 
corrections for the reactor angle of the order 613 ~ A^, where Ac ~ 0.2. In this paper, we discuss 
the possibility to achieve large 6113 within the T2K region with maximal atmospheric mixing angle, 
sin 623 ~ 1/2, and trimaximal solar mixing angle, sin 6\_2 = 1/3, through the deviation from the 
exact tri-bimaximal mixing. We derive the structure of neutrino mass matrix that leads to the large 
6\z leaving maximal 62Z and trimaximal 612- It is shown that such a structure of neutrino mass 
matrix can arise in a model with Si flavor symmetry. 
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I. INTRODUCTION 

(~| . The T2K collaboration recently gives that the reactor angle is [l| 

0.087(0.100) < sin6ii3 < 0.275(0.306), (1) 






with best fit value of s\n6i2 = 0.17(0.19) for normal (inverted) hierarchy in the neutrino masses. Such a result must 
be taken as a hint since it is only at 2a. We however consider the implications of such an important indication in this 



■ paper. 
\jj ■ In the last decade, the tri-bimaximal (TB) mixing pattern introduced by Harrison et al. in 2002 [2| 

o 



V2 



has been used as a guide in neutrino physics for the fiavor problem. However if the result of T2K will be confirmed, 

I ^ it will have strong impact on this point of view. In fact, most of the models predicting tri-bimaximal mixing at the 

:• ' leading order are compatible only with small values of the reactor angle. In a generic model, the three mixing angles 

.5^ ' receive corrections of the same order. Departure of the solar angle from the trimaximal values are at most of ©(A^) 

S^ . where Ac ~ 0.2. Therefore it is possible to have a deviation of order A^ « 0.04 only Q for the reactor angle which 

$^ ■ is about half of the lower bound in ([l} . This is only an estimation and it must be considered individually case by 

. . 1 case. We however expect that most of them are on the border of validity if not excluded completely. Therefore it 

is important to search for the models with large ^13 and with maximal atmospheric mixing angle and trimaximal^ 

solar angle. Recently, some works fitting T2K result has been presented Q. There are some models based on discrete 

fiavor symmetries before T2K data which predicts large reactor mixing angle, for an incomplete list see reference [y], 

and for a classification of models with fiavor symmetries classified by its predictions for reactor angle see [7| 

In this paper, we study the possibility to obtain the large reactor angle with tri-bimaximal values of the solar and 
atmospheric mixing angles. The lepton mixing matrix with such mixing pattern was first proposed by King in [8| 
and called Tri-bimaximal-reactor (TBR) mixing. Using such mixing matrix, we found the structure of the deviations 
in the neutrino mass matrix from its TB texture which leads to TBR mixing. We then show that such a particular 
deviation in neutrino mass matrix can arise in a model with Si flavor symmetry. 
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^ The word "trimaximal" is used for difi'erent meaning in previous studies (jj. We mean here sin^ 612 = 1/3 by trimaximal solar angle. 



The paper is organized as follows. We discuss the conditions to obtain a mass matrix with maximal atmospheric 
mixing angle, trimaximal solar mixing angle and a non-zero reactor mixing angle within the T2K region in section 
HTl In section ITnl we present a model based on the group of permutation of four objects, S4 where the neutrino mass 
matrix with particular form discussed in section |TT] is obtained. Finally, we discuss the phenomenology of the model 
and conclude in in section IIVI 



II. LARGE REACTOR TRI-BIMAXIMAL MIXING AND NEUTRINO MASS MATRIX 

In this section, we study the structure of the neutrino mass matrix (in the diagonal basis of the charged leptons) 
that gives maximal atmospheric angle ^23 = 7''/4, trimaximal solar angle sin0i2 = l/vS and an arbitrary reactor 
angle 613 = A. In the standard PDG [9| parametrization, the lepton mixing matrix with the above values of mixing 
angles is given by [8] 



TT — 

UtBR = ^23( j) ^13(A) -^12(^12) = 
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0{X'). (3) 



We do not consider the CP violation in the lepton sector assume that the above parameters are real for simplicity. 
The neutrino mass matrix diagonalized by ([3]) is given by 



m 



Utbr ■ mt'^s ■ f/^BR = rnl"" + 5m, (4) 



where mf,^^^ is a diagonal matrix with the neutrino mass eigenvalues, m,^, m,^ and m,^. This leads to the following 
structure of the neutrino mass matrix 

2y — X X X 

rn^^ = \ X y + z y- z \ , (5) 

X y — z y + z 

where x — (to2 — "t-i)/3, y = (mi + 2m2)/6 and z ~ m^/2 and 



dm, = A ai /3i + aM a2 /32 -/32 + y A" a„ , (6) 
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with ai = — (x— 2y + 2z)/\/2, /3i = \/2x, 0.2 — —x/2, (i2 — —{x— 2y + z)/2 and 7 = 2; — 2y + 2z. Note that /32 can be 
reabsorbed into the TB term m^^. The above form of neutrino mass matrix predicts maximal atmospheric mixing 
angle and trimaximal solar mixing angle if all the terms with all powers of A are taken into account. If one truncates 
the series in eq. © at n < 3, the neutrino mass matrix then implies 

• (A) negligible deviations from maximality in the atmospheric mixing angle; 

• (B) small deviation from trimaximality in the solar mixing angle; 



• 



(C) prediction of Oi//3/3 ex A 



2 



The prediction (C) is evident from eq. ([6]) and we verify numerically (A) and (B) in the section HVl We observe that 
the main structure of the deviation Sm, of order A in eq. ^ is ^-t antisymmetric, see [lC|^. Therefore a possible 
flavor symmetry with neutrino mass matrix texture ^ must contain the group S2 of the fi-r permutation and must 
be compatible with tri-bimaximal in the unperturbed limit. One possible flavor symmetry with such features is S4 
which contains 5*2 as a subgroup and leads to tri-bimaximal mixing [ll| . 



^ Note that the main structure of the deviation 5m^ of order A in eq. JgJ is similar to the one found in the paper by T. Araki in Ref. [3I 
where (contrary with respect to us) the solar angle is not fixed to be the trimaximal one. 



III. THE MODEL 

We assume 5*4 (see appendix) flavor symmetry and extra abelian Zjsi symmetry in order to separate the charged 
leptons from the neutrino sector as usual in models for TB mixing, see for instance [3|. In order to simplify the model 
as much as possible and to render more clear the main features of the model, we do not enter into the details of the 
particular Zjv symmetry required in this model. Our purpose is to show that the neutrino mass matrix ((4]) with the 
structure given by (O and ^ can be obtained from symmetry principle. We assume that light neutrino masses arise 
from both type-I and type-II seesaw and introduce only one right-handed neutrino. The matter content of our model 
is given in tableU) 
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TABLE L Matter content of the model giving TB mixing at the leading order 



In the scalar sector, we have one SUl{2) triplet A and one singlet (j) in the neutrino sector transforming both as 3i 
of 5*4. We have two electroweak singlets (pi and £^1 in the charged lepton sector, transforming as doublet and singlet 
of Si respectively. As it has been already mentioned, the two sectors can be separated by introducing an abelian Zn 
symmetry under which /^, ipi and £1 are charged while the other fields could be singlets of Zjsf. The Yukawa interaction 
of the model is 



'^i = -ryi{LlR)i,h£,i + —y2{LlB,)2lT''Pi + h.c. 



(7) 
(8) 



where A is an effective scale. We assume the following S4 alignment in the vacuum expectation values (vevs) of the 
scalar fields. 



(A")-i;a(1,1,1)^, (0)=«^(O,l,-lf, {^)^{v,,V2f 



(9) 



where vi 7^ V2 ■ Using the product rules shown in appendix A, one can easily see that the charged lepton mass matrix 
is diagonal and the lepton masses can be fitted in terms of three free parameters yi, vi and V2, see |12| for details. 

The type-II seesaw gives a contribution to the neutrino mass matrix with zero diagonal entries and equal off diagonal 
entries since it arises from the product of three S'4 triplets. Since we introduced only one right-handed neutrino, Dirac 
neutrino mass matrix is a column m/) ^ (0, 1, —1)"^ and the light-neutrino mass matrix from seesaw relation is given 
by 



m'yP^-'^-mDml 




(10) 



Considering both the type-I and type-II contributions, we have the light neutrino mass matrix which can be diago- 
nalized by TB mixing matrix |13|^ 



m 



TB 




(11) 



The mass eigenvalues of the above matrix are mi — —a, 7712 = 2a and 7713 = —a + 26. Here a — yaVA and 
b — j/^w^w?/(A^M) where Vh — (^°)- This neutrino mass matrix is compatible with the normal hierarchy only and 
predicts zero neutrinoless double beta decay rriee = 0. 



^ In [13|| . similar structure l llOII has been obtained through only type-II seesaw and S4, symmetry. 



In order to reproduce deviations like eq. ([6]) in the neutrino mass matrix, we introduce in the scalar sector one Higgs 
triplet Arf that transforms as a doublet under 5*4 and an electroweak singlet (jj^ that transforms as a triplet 3i under 
S4. With inclusion of these fields, the Yukawa interaction Lagrangian Ci, contains also the terms 

- £, D VfiLLAd + ^{L(i3d)uhvR + h.c. (12) 

We assume that A^ and (t>d take vevs along the following directions 

(A°)=i;d(l,Of, ((/.rf)=wd(l,0,0f. (13) 

Here we also assume that ya.(} <C Ha.b- This can be realized assuming that A^; and (pd are charged under some extra 
abelian symmetry like Z^r or Ufn{^)- 

After electroweak symmetry breaking and integrating out the right-handed neutrino, eq. (J12l) gives the following 
contribution to the neutrino mass matrix 

^g^(^^)l,(^</)d)i, + |||:(^0d)l,(z/0d)l,. (14) 

The second term in eq. (1141) is smaller with respect to the first since we have assumed j/q, -^ yt- In particular assuming 
?/b ~ 1 and ya ^ ^ the first term is proportional to A and the second term is proportional to A^ . The extra contributions 
to the neutrino mass matrix from the type-I see-saw are as follows 

Sm'JP°-'^ci\\ 10 +C2A2 I I . (15) 
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where, ci and C2 are coefficients of order 0(1). From the extra type-II seesaw term in eq. (fT2|) and using the vev 
alignments as in (jlSp. the additional contribution to the perturbed neutrino mass matrix will be proportional to 
vii^i — V2Vi^ therefore the contribution to the neutrino mass matrix coming from Type-II see-saw is 




Jm*/P^-" - 1 . (16) 

Vo -x) 

Putting all these results together, the structure of the deviation in neutrino mass matrix can be written is 



dull, 



( i a' 


-a' 


a' /3' 
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(17) 



where a' = ybVavlv^Ud/iA'^M), fi' = ypvd, 7' = vlvlv^Ud/iA'^M). 

The deviation obtained in our model equal to the neutrino mass deviation in eq. Qi truncated at A^ with a2 = 0. 
Such a difference does not modify significantly the prediction of maximal atmospheric angle and trimaximal solar 
angle. In the next section, we study the phenomcnological implication of our neutrino mass texture. 

IV. PHENOMENOLOGY 

Combining eq. ()lip and eq. (|17|) . the resulting neutrino mass matrix in our model is 

(18) 

As mentioned earlier, we assume ya,yb ^ 0(1) and ya,y/3 ^ 0(A) which implies the hierarchies in the elements 
a,b^ a',P' :S> 7' in the above structure. Since the {m^)ii entry is 7 ^ O(A^), we have neutrino less double beta 
decay rate iriee oc A^ and for small values of A as in the unperturbed case only the normal neutrino mass hierarchy can 
be fitted. The neutrino mass matrix (jl8|) obtained from the model is equivalent to the matrix ^ up to the correction 
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of ©(A^), the neutrino mass matrix p8|) is diagonalized by the mixing matrix ([3]) up to 0{X^) corrections. Note that 



the neutrino mass matrix (J18p obtained in the model has 5 free parameters while the derived structure in eq. ([3]) has 4 
real parameters (x, y, z and A). We fix free parameters of eq. (|18|) in terms of the parameters of eq. ([3|) by comparing 
both the structures at each order of A. Comparing the leading order expressions of the neutrino mass matrix in 
eq. (ITTt and eq. (O, we restrict our parameters to be 



X = 2y; a = 2y; b^y + z. 
Further comparing the higher order terms in A, we obtain the relations 

a' = -V2zX; /?' = 2V2yX; 7' = 2zX^ 



(19) 



(20) 



Note that this is not the most general case for our model, nevertheless we want to point out that even in the case of 
bimaximal atmospheric mixing angle and trimaximal solar mixing angle it is possible to obtain large reactor mixing 
angle. Also, we expect the negligible deviations from the tri-bimaximal values for the solar and atmospheric mixing 
angles due to the fact that our model predicts a2 = if compared with eq. Q and it generates the terms only up to 
the 0{X^). We analyze such deviations by randomly varying y, z and A with the constraints that the square mass 
differences and the mixing angles are in the observed (Scr) range of validity [IJ]. The results of our analysis are shown 
in figure [T] It is evident from figure [1] that for restricted parameter space as specified earlier, model allows large 
^13 with negligible deviations in atmospheric and solar mixing angles from there bimaximal and trimaximal values 
respectively. We also check the predictions for the neutrinoless double beta decay rates in our model for restricted 
parameter space specified above and find that the region for 4.5 meV < m^i < 5.8 meV and 0.5 meV < \mpj3\ < 3.5 
meV is allowed for the values of ^13 in the 2a limits indicated by T2K. 
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FIG. 1: The figure in the left side shows the allowed region for ^13 vs. 623- The figure in the right side shows the allowed region 
for 9is vs. 612. The horizontal continuous (dashed) hues represent the best fit value {2a deviations) of T2K for the reactor 
neutrino mixing angle. 



In summary, we found the structure for the deviation in the neutrino mass matrix from the well known TB pattern 
in such a way that the lepton mixing matrix has large atmospheric mixing angle and trimaximal solar mixing angle 
with an arbitrary large reactor angle. The deviation must be approximately /z-r antisymmetric. This fact suggests 
us that the flavor symmetry could be some permutation symmetry containing S2 {fJ^-T exchange) subgroup. S3 is too 
small since it does not give the TB mixing. The smallest permutation group with this property is S4. We provide a 
candidate model based on 5*4 where in the unperturbed limit the neutrino mass matrix is TB. Then assuming extra 
scalar fields we show the possibility to generate deviations from the TB that give a large ^13 in agreement with T2K 
result, maximal atmospheric mixing angle and trimaximal solar mixing angle in good agreement with neutrino data. 
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Appendix A: 54 product rules 

In the basis where the generator of 5*4 are real, the products oi fj, x ^ (see [lll|): 
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and for 2 x 82 
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